setwd("D:/project/GGT210524001RS-24")


info <- read.delim(file="vali.txt",header = T,row.names = 1,sep="\t")
info$lghr <- log(info$hr)
lglci <- log(info$lci)
lguci <- log(info$uci)
info$selghr <- (lguci-lglci)/(2*1.96)

library(meta)
os <- metagen(TE = lghr,seTE = selghr,studlab = rownames(info),data = info,sm = 'HR',comb.fixed=FALSE)
pdf("meta.pdf",height=3,width=8)
forest(os,layout = "revman5")
dev.off()
forest(metainf(os,pooled = 'random'),layout = "revman5")


# vali.txt
# 	pvalue  hr	lci	uci
# GSE10886	0.999	19	0	99999
# GSE10893	0.223	0.402	0.118	1.374
# GSE45725	0.384	2.733	0.409	18.27
# GSE53031	0.298	0.479	0.150	1.533
